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PROCEDE POUR I* 1 ESTIMATION 
DES C ARACTERI ST I QUE S D'UNE PRECIPITATION 



La presente invention concerne un procede pour 
5 1 1 estimation des caracteristiques d'une precipitation, et en 
particulier du taux precipitant pour une precipitation solide. 

On connait dans 1 1 etat de la technique 1 1 utilisation des 
caracteristiques d'une image radar pour estimer les 
caracteristiques pluviometriques . En particulier, le brevet 
10 europeen EP1049944 decrit une technique pour 1 1 estimation de 
pluie avec un radar. Selon ce brevet, on realise les etapes 
suivantes : 

-on mesure, au moyen dudit radar bipolaire, sur un 
intervalle [ro, ri] donne de rayon de trajet par 
15 rapport audit radar, la phase dif f §rentielle (Odp et 

la reflectivite apparente Z. selon au moins l'une des 
polarisations H ou V ; 

- on determine une estimation de la valeur No* 
representative de la distribution dimensionnelle des 

20 gouttes de pluie, a partir de la difference de . la 

phase dif f erentielle entre ro et ri et a partir d'une 
integrale d'une fonction de la reflectivite apparente 
Z, le long de 1 1 intervalle [ro, ri] ; 

- on deduit la valeur du taux de precipitation en un 
25 point a partir de No* et de la reflectivite apparente 

en ce point. 

On connait egalement un precede pour 1' estimation de 
pluie decrit dans le brevet PCT WO03007016 decrivant un 
procede pour l f estimation d'un taux precipitant au moyen d'un 
30 radar '.bipolaire, caracteris<§ par les diff^rentes etapes 
suivantes : 
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- on mesure au moyen dudit radar bipolaire, sur un 
intervalle [rl, rO] donne de rayon r de trajet par 
rapport audit radar, la phase dif f erentielle Odp et 
la reflectivite attenuee Z selon au moins l'une des 
polarisations H ou V ; 

-on determine une estimation de la valeur K(rO) de 
1 ' attenuation specif igue en rO a partir du profil de 
reflectivite attenuee ainsi mesure , ainsi qu'a partir 
de la difference de la phase dif f erentielle entre rO 
et rl; on determine une estimation K(r) de 
1' attenuation specif ique en r en fonction de 
1 ' attenuation K(rO) ainsi determined et du profil de 
reflectivite attenuee Z(r); on determine le taux de 
precipitant R(r) connaissant K(r) . 
Ces differentes solutions permettent de caracteriser des 

precipitations liquides et d'estimer de fagon precise le taux 

de pluie (en mm/h) , mais pas d'estimer les caracteristiques 

des precipitations solides telles que la neige. 

L' invention concerne un nouveau procede visant a remedier 

a cet inconvenient. 

A cet effet, 1 ' invention concerne selon son acception la 
Plus generale un procede pour 1' estimation des 
caracteristiques d ' une precipitation, et en particulier du 
taux precipitant pour une precipitation solide comportant une 
etape d ' acquisition d 1 une image radar comportant au moins un 
plan vertical d'une zone de precipitation et d'un traitement 
d'un profil vertical pour fournir des signaux numeriques 
representatifs de la reflectivite selon la direction verticale 
h caracterise en ce que 1'on realise une etape d ' integration 
desdits signaux representatifs de la reflectivite pour 
delivrer un signal representatif du profil dans le plan 
vertical du diametre moyen des particules pondere par la masse 
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de chaque particule, et une etape de determination de la 
concentration des particules solides a partir des signaux 
calcules dans les etapes precedentes . 

De preference, 1'etape d ' integration consiste a 
5 determiner la variable Z(h) l'observable radar en mm 6 /m 3 en 
fonction de 1 ' altitude h a partir de ladite image radar, et a 
determiner ledit diametre moyen des particules Dm(h) par 
resolution de 1' equation : 

J^=-0.25^ # aD^ 5 10- l8 Z+|^^)). (2) 

10 ou : 

■ Z est l'observable radar a inverser en mm 6 ra" 3 ; 

■ D m est en m ; 

■ a et b sont des coefficients specif iques des 
particules de type « agregats » . Le coefficient a est par 

15 exemple egal a 35184 et le coefficient b est egal a 3.16- 

■ k Gff est le coefficient d'efficacite du processus 
d'agr^gation a ajuster, ce coefficient k eff etant par 
exemple egal a 0*3 

L ' integration de (2) necessite une condition a la ldLmite 
20 d 1 integration • Avant ageu semen t, ladite condition a la limite 
d' integration est determinee afin que la valeur D m (h) au sommet 
du nuage corresponde a la valeur predeterminee du nombre total 
de particules au sommet du nuage. 

Selon un mode de realisation avantageux, on determine le 
25 profil du nombre total de particules n t (h) par l'equation 
suivante : 

n T (h)=x.Z(h)/D ra (h) 6 

ou x est egal a 25,4 lO^ 18 

Selon un autre mode de realisation, on determine le 
30 parametre meteorologique N 0 (h) par l'equation suivante : 
N 0 (h)=y.Z(h)/D m (h) 7 
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Ou y est egal a 102 10" 
Selon un troisieme mode de raise en oeuvre, on determine 
le parametre meteorologique correspondant au profil du contenu 
en glace IWC(h) [en g/m 3 ] par 1' equation suivante : 
IWC(h)=wZ(h)/D ra (h) 3 
ou w est egal a 1,25 10" 12 

Selon un quatrieme mode de mise en oeuvre, on determine 
le parametre meteorologique correspondant au profil du taux 
precipitant R(h) solide (mm/h equivalent fondu) par 1' equation 
suivante : 

R(h)=r.Z(h)/D m (h) 2 - 35 

Ou r est egal a 4,698 10' 10 

L' invention sera mieux comprise a la lecture de la 
description qui suit, se referant a un exemple non limitatif 
de realisation. 

Le precede s' applique essentiellement aux precipitations 
stratif ormes . il considere que les noyaux glagogenes ne sont 
actives qu'a temperature fortement negative, c'est-a-dire au 
sommet du nuage. Les cristaux de glace formes a haute altitude 
sedimentent et grossissent au cours de leur chute soit par 
sublimation de la vapeur d'eau saturante ambiante, soit par 
collection et givrage de gouttelettes d'eau nuageuse 
surfondues, soit par agregation au hasard de leurs collisions 
avec d'autres cristaux de glace. Des trois processus de 
croissance, seule l'agregation change la concentration en 
particules de glace. Le precede de caracterisation repose . 
essentiellement sur une description simplifiee du mecanisme 
d' agregation. L ' invention est basee sur une classe de procedes 
« profileurs », en ce sens qu'elle inverse le profil vertical 
de reflectivite mesure dans la glace, pour en tirer le profil 
vertical du taux precipitant solide. 

Les etapes de la methode d' inversion sont les suivantes : 



1- La distribution dimensionnelle des particules 
exprimee en « diametre equivalent fondu » , est supposee 

exponentielle i.e. : 

N(D)=N 0 exp(-4D/Dm) (1) 
5 ou N(D) est la concentration en particules par m 3 et 

par intervalle de diametre , et N 0 et D m sont les deux 
parametres qui caracterisent la distribution. 

2- On determine le sommet h max et la base h min de la 
couche de precipitation solide. 

10 a- h max est 1' altitude maximum du profil de 

reflect ivite mesure Z(h) . 

b. est soit l f altitude de 1'isotherme 0°C si 

la temperature au sol est positive, soit le niveau du 

sol si la temperature au sol est negative, 
15 3- On determine alors le profil du parametre D m (h) 

entre et h ^ en resolvant 1' equation 

dif ferentielle : 



f^^0.25k e/r aDt 5 lO^Z ^ | £| }) m ( 2 



) 



20 ou 



Z est 1 'observable radar a inverser en 



mmV 3 ; 

D m est en m ; 

a et b sont des coefficients specifiques 
25 des particules de type « agregats », egaux 

respectivement a 35184 et 3.16 a partir des 
observations de Locatelli et Hobbs (1974) ; 

k e£f est le coefficient d'efficacite du 
processus d'agr^gation a ajuster (la valeur k eff 
30 =0.3 semble correcte ) . 
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4- L' integration de (2) s'effectue a partir du 
haut, ou la condition a la limite s'exprime en fixant le 
nombre total de particules n T (ou nombre de noyaux 
glagogenes actives au sommet de nuage). On peut prendre 
n T (h aax )~ 10 6 m~ 3 , ce qui permet d'exprimer la condition a 
la limite D m (h max ) comme : 

Dm(H max )=25,4 10r' a (Z(h max /nj{h mc J) I/6 ( 3 ) 

5- Le profil D m (h) entre h aax et h ^ une fois 
determine, on calcule les profils des autres parametres 
d'interet par les expressions suivantes : 

a. Profil de N 0 : N 0 (h)=102 . 10" I2 Z (h) /D m (h) 7 

b. Profil du nombre total de particules n T (h) 
[en m" 3 ] : N T (h)=102.10- 12 Z(h)/D ro (h) 6 

c. Profil du contenu en glace IWC(h) [en g/m 3 ] 
: IWC(h)=l, 25.10- 12 Z(h)/D ra (h) 3 

d. Profil du taux precipitant R(h) solide 
(mm/h equivalent fondu) 

En utilisant la loi de vitesse terminale de 
chute determinee par Locatelli and Hobbs pour des 
agregats: [v r =107.6 D 0 65 (D in m) ] , R(h) s'exprime 
par : 

R(h)=4 / 698.10- 10 Z(h)/D m (h) 2 ' 35 
La figure 1 represente un exemple de profil vertical de Z 
a inverser (dans cet exemple, l'isotherme 0°C est au niveau du 
sol) . 

La figure 2 represente la comparaison du profil de D m 
resultant de 1' inversion de Z par le modele d'agregation avec 
l'estimateur classique. 

La figure 3 represente la comparaison des profils de N 0 et 
n T resultant de 1' inversion de Z par le modele d'agregation 
avec l'hypothese classique et des observations. 
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La figure 4 represente la comparaison du profil de R 
resultant de 1' inversion de Z par le modele d'agregation avec 
l'estimateur classique. 



4 
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REVENDICATIONS 



1 - Procede pour 1 ' estimation des caracteristiques d'une 
precipitation, et en particulier du taux precipitant pour une 
precipitation solide comportant une etape d ' acquisition d'une 
image radar comportant au moins un plan vertical d'une zone de 
precipitation et d ' un traitement d'un profil vertical pour 
fournir des signaux numeriques representatif s de la 
reflectivite selon la direction verticale z caracterise en ce 
que l'on realise une etape d ' integration desdits signaux 
representatifs de la reflectivite pour delivrer un signal 
representatif du profil dans le plan vertical du diametre 
moyen des particules pondere par la masse de chaque particule, 
et une etape de determination de la concentration des 
particules solides a partir des signaux calcules dans les 
etapes precedentes. 

2 - Procede selon la revendication 1, caracterise en ce 
que 1' etape d ' integration consiste a determiner la variable 
Z(h) 1' observable radar en mmVm 3 en fonction de 1' altitude h a 
partir de ladite image radar, et a determiner ledit diametre 
moyen des particules Dm(h) par resolution de 1 'equation : 

(2) 



OU 



Z est 1' observable radar a inverser en 

mm 6 m" 3 ; 

■ D m est en m ; 

a et b sont des coefficients 
specifiques des particules de type 
« agregats » . 
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■ k e£f est le coefficient d'efficacite 
du processus d'agregation a ajuster. 

3 - Procede selon la revendication 2, caracterise en ce 
5 que le coefficient k eff est egal a 0.3. 

4 — Procede selon la revendication 2 ou 3, caracterise en 
ce que le coefficient a est egal a 35184. 

10 5 — Procede selon la revendication 2, 3 ou 4, caracterise 

en ce que le coefficient b est egal a 3.16. 

6 — Procede selon 1'une quelconque des revendications 2 a 
4, caracterise en ce que ladite constante d' integration est 

15 determinee afin que la valeur D m (h) au sommet du nuage 
corresponde a la valeur predeterminee du nombre total de 
particules au sommet du nuage. 

7 — Proc6d6 selon l r une quelconque des revendications 
20 precedentes , caracterise en ce que l'on determine le profil du 

nombre total de particules n^h) par 1' equation suivante 'z 
n T (h)=x.Z(h)/D ro (h) $ . 

8 — Procedg selon la revendication precedente, 
25 caracterise en ce que x est egal a 25,4 10" 18 . 

9 — Procede selon l'une quelconque des revendications 
precedentes, caracterise en ce que l f on determine le parametre 
meteorologique N 0 (h) par 1' equation suivante : 

30 N 0 (h)=y.Z(h)/D m (h) 7 . 
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10 — Procede selon la revendication precedente, 
caracterise en ce que y est egal a 102 10" 1 
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11 - Procede selon l'une guelconque des revendications 
precedentes, caracterise en ce que l'on determine le parametre 
meteorologique correspondant au profil du contenu en glace 
IWC(h) [en g/m 3 J par 1' equation suivante : 

IWC(h)=wZ(h)/D m (h)\ 

12 - Procede selon la revendication precedente, 
caracterise en ce que w est egal a 1,25 10" 12 . 

13 - Procede selon l'une quelconque des revendications 
precedentes, caracterise en ce que l'on determine le parametre 
meteorologique correspondant au profil du taux precipitant 
R{h) solide (mm/h equivalent fondu) par l'equation suivante : 

R(h)=r.Z(h)/D m (h) 2 ' 35 . 

14 - Procede selon la revendication precedente, 
caracterise en ce que r est egal a 4,698 10 -10 . 
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